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The cross-correlation between the integrated Sachs- Wolfe (ISW) effect and the large scale structure 
(LSS) is a powerful tool to constrain dark energy and alternative theories of gravity. In this paper, 
we obtain observational constraints on kinetic gravity braiding from the ISW-LSS cross-correlation. 
We find that the late-time ISW effect in the kinetic gravity braiding model anti-correlates with large 
scale structures in a wide range of parameters, which clearly demonstrates how one can distinguish 
modified gravity theories from the ACDM model using the ISW efi'ect. In addition to the analysis 
based on a concrete model, we investigate a future prospect of the ISW-LSS cross-correlation by 
using a phenomenological parametrization of modified gravity models. 
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I. INTRODUCTION 



Modifying general relativity at long distances is a pos- 
sible way to explain the present accelerated expansion of 
the Universe indicated by different cosmological observa- 
tions such as type la supernovae [l|, 0] , the cosmic mi- 
crowave background (CMB) anisotropies the large 
scale structure (LSS) of galaxies 0,0, and clusters of 
galaxies 0, H- ^ number of modified gravity models 
have been proposed so far, including scalar-tensor the- 
ories, f{R) gravity I9l4l5i. the Dvali-Gabadadzc-Porrati 



(DGP) brane model [16|, |17[, and galileon gravity jl8j . 

Recently, the galileon gravity model has attracted con- 
siderable attention, which is constructed by introduc- 
ing the scalar field with the self-interaction whose La- 
grangian is invariant in the Minkowski space-time under 



the Galilean symmetry 9^ 



dnfj) + bn, which keeps 



equation of motion the second order differential equation 
with no ghostlike instabilities [l^- A remarkable prop- 
erty of the galileon model is the Vainshtein mechanism 
that the self-interaction terms induce the decoupling of 
the galileon field (j) from gravity at small scale [T^ . This 
allows the galileon theory to recover general relativity 
around a high density region, which ensures the consis- 
tency with the solar system experiments. Furthermore, 
the galileon model can lead to the late-time accelerated 
expansion of the universe, which can be useful as an al- 
ternative to the simple dark energy model. 

Cosmology of the galileon model has been extensively 
investi gate d by many authors pol - lssj . Recently, Deffayet 
et al. [33 proposed the most generalized scalar-tensor 
theory, including f{R) theories, galileon theories and the 
kinetic gravity braiding j40l-l42j . In addition, several au- 
thors have investigated the observational constraints on 
various galileon theories from type la supernovae, the 
CMB shift parameter, baryon acoustic oscillations and 
the growth rate of matter density perturbations [43l - l46j . 

It was pointed out that the integrated Sachs- Wolfe ef- 
fect can distinguish between the galileon model and the 



ACDM model, which was discussed in the scalar-tensor 
galileon model [23|. In contrast to the ACDM model, in 
which the gravitational potentials decay because of the 
accelerated expansion of the universe, the amplitude of 
the gravitational potentials in the galileon model may in- 
crease due to perturbations of the galileon field. In these 
two models, the integrated Sachs- Wolfe effect predicts 
opposite signs. Therefore, this provides us with a useful 
chance to test the galileon model using a cross-correlation 
between the matter (galaxy) distribution and the cos- 
mic microwave background anisotropies \at\ . l48j . In the 
present paper, we focus our investigation on the observa- 
tional constraints from the cross-correlation between the 
CMB anisotropies and the galaxy distributions. 

Throughout the paper, we use units in which the speed 
of light and the Planck constant are unity, c = h = 1, and 
Mpi is the reduced Planck mass related with Newton's 
gravitational constant by Mp\ = 1/ vSttG. Wc follow the 
metric signature convention (— , +, +, +). 



II. KINETIC GRAVITY BRAIDING MODEL 

We start with the most general minimally coupled 
scalar field theory with second-order field equations [Sll, 
|40||. whose action is of the form 



5" = 



A/2 



R + K{4>, X) - G{4>, X)n(j> + 



,(1) 



where R is the Ricci scalar, K{(I>,X) and G{4),X) are 
arbitrary functions of the scalar field and its kinetic 
term X := -.g^''V,,0V^0/2, Ucf) = ^^^^V^V^^, and 
is the matter Lagrangian. Although is not directly 
coupled to R, the interaction GDc/) causes "braiding" of 
the scalar field and the metric, giving rise to non-trivial 
and interesting cosmology. In this paper, we consider 
the shift symmetric case, K = K{X) and G = G(X), 
for which an attractor solution is generically present [3J, 
|40| . leading to a late-time de Sitter expansion driven by 



2 



constant kinetic energy 
the following model [4^ 

KiX) 



of (f). In particular, we focus on 



G{X) = Ah 



PI 



X 



(2) 
(3) 



where n and Vc are the model parameters, though Tc is 
required to be tuned as Tc ^ ^^(7^- Throughout the pa- 
per we consider for simplicity the cosmological evolution 
along the attractor, assuming that the background solu- 
tion converges to the attractor at some sufficiently early 
time. 

We work in the spatially flat Friedmann-Robertson- 
Walkcr Universe and consider the metric perturbations 
in the conformal Newtonian gauge. 



,(4) 



where the scale factor a is normalized to unity at the 
present epoch, rj is the conformal time, x is the co- 
moving distance and c?f2(2) the line element on a unit 
sphere. The quasi-static approximation is applicable for 
0{Pcl/o?) > 0{H'^), which corresponds to n < 10 for 
the wavenumber k ^ 0.01/iMpc^^. Under this approxi- 
mation, the modified Poisson equation may be written 



(5) 



where pm is the energy density of matter, 6 is the matter 
density contrast, and Gcff is the effective gravitational 
coupling, 



GcS = G 



2n -\- 3nrii„(a) — fi,„(a) 



(6) 



f^m(a)[57i — i7in(a)] 

The traceless part of the Einstein equations is given by 

* + $ = 0. (7) 

The evolution equation for the matter overdensity 5 is 
given by 



5 + 2HS : 



(8) 



For large n, the quasi-static approximation no longer 
works. In this case, one needs to solve the full pertur- 
bation equations. The full perturbation equations and 
the validit y o f the quasi-static approximation are inves- 
tigated in [42l. For convenience, we define Ukir]) as 



- k 



2 *k - $k 3 



(9) 



where VPk, 'fk and (5k are the metric perturbations and 
the energy density contrast in Fourier space, respec- 
tively. As long as the quasi-static approximation holds, 
Ukiv) = (Geff/G)(i?i(r?)/ai^i(r?o)), where is the 



growth factor of the kinetic gravity braiding model nor- 
malized as Di(rf) = a at early stage of the evolution and 
770 is the present conformal time. 

The authors of Ref. [l^ obtained the observational 
constraints on the kinetic gravity braiding model with 
the functions ([2]) and ([3]) from type la supernovae, the 
WMAP CMB anisotropy experiment, and the Sloan Dig- 
ital Sky Survey (SDSS) luminous red galaxy (LRG) sam- 
ples. The result of Ref. [12] indicates that models with 
larger n better fit the CMB distance observation, though 
the constraints derived there are not so strong. In this 
paper, we will improve the constraints significantly by 
using the galaxy-CMB cross-correlation. 



III. ISW-GALAXY CROSS-CORRELATION 

The ISW effect contributes to the power spectrum of 
the CMB anisotropy on large scales. However, the direct 
detection of the ISW effect is very difficult due to the 
cosmic variance and the relatively small amplitude of the 
ISW effect in the power spectrum. Nevertheless, by cross- 
correlating the CMB anisotropics and the galaxy number 
density fluctuations, one can isolate the ISW effect, which 
will be a powerful probe of dark energy and modified 
gravity models. The ISW effect in the CMB anisotropy 
is given by 



Ar(7) 

T 



dr;[v?'(r;,x)-$'(ry,x)]. 



(10) 



iid 



where rjd is the decoupling time, which may be taken to 
be zero practically, and the prime denotes the differenti- 
ation with respect to r/. The fluctuations in the angular 
distribution of galaxies is given by 



N 



dzdg{Tj,:K.)Wiz), 



(11) 



where yV(x) is the selection function and Sg is the galaxy 
number density contrast, which is related to the energy 
density contrast S through the bias b as Sg = b{k, z)S. 
The cross-correlation of the galaxy fluctuations and the 
CMB anisotropics is thus expressed as 



AT(7) A7Vg(f ) 



T 



N 



1 

in 



J2i'2l+l)CiVi{fi), (12) 



where Ve is the Legendre polynomial with fj. being the 
cosine of the angle between 7 and 7' , 



(^+ 1/2)2 
dUkiv) 



dzH{z)W{2 



Di{z = 0) 



dz 



b{z,k)P(k) 



(13) 



k={l+l/2)/x 



and P{k) is the matter power spectrum at the present 
time, rj = rjQ. In deriving the above equation, we used 
the small angle approximation, \. 
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IV. COMPARISON WITH OBSERVATIONS 

Figure [T] shows the measured cross-correlation func- 
tion for the six catalogues (2MASS, SDSS galaxies, LRG, 
NVSS, HEAO, and QSO) available in Ref. [H. ^ The 
solid and dashed curves show the theoretical predictions 
for the ACDM model and the kinetic braiding model with 
n = l, 10, 100, 1000, and 5000. We fixed the baryon den- 
sity as fib = 0.0451, the matter density as floh"^ = 0.1338, 
the bubble constant as /i = 0.702, the spectral index as 
Us = 0.966, and the amplitude of the density fluctuation 
asA^ = 2.42xl0~3at A:o = 0.002 Mpc"! [1]. The galaxy 
bias of each catalog is found in Ref. [5^ and the values 
are 6 = 1.4, 1.0, 1.8, 1.5, 1.06 and 2.3 for 2MASS, SDSS 
galaxies, SDSS LRG, NVSS, HEAO, and QSO, respec- 
tively. For the theoretical modeling of the power spec- 
trum, we used the linear power spectrum with the trans- 
fer function T{k) of Eisenstein and Hu |76[ and obtained 
$, and S by solving the full perturbation equations. 
As we can see from Fig. [1] the cross-correlation function 
for the kinetic gravity braiding model with larger n is 
closer to that for the ACDM model, and they are practi- 
cally indistinguishable for n > 0(1000). This is because 
the perturbation of the kinetic term 5X become zero in 
the limit of ti — > cx) [4^ . Another important feature 
found in Fig.[T]is the anti-correlation for the kinetic grav- 
ity braiding model with small n. This is caused by the 
enhancement of the effective gravitational coupling GcS 
and the consequent growth of the gravitational potential, 
leading to the opposite sign of the function dUk{z)/dz 
compared to the ACDM model. Since the measurements 
of the cross-correlation function in each catalog shows a 
positive correlation, kinetic gravity braiding with small 
n will be inconsistent with observations. 

In our analysis, the total chi squared is given by 



Xtotal 



,obs 



„thco\ 



)C,^i(cf^-cf™), (14) 



where c°^^ is the cross-correlation function obtained from 



„thco 



oretically predicted from Eq. p^ . and Q-^ is the in- 



verse of the covariance matrix obtained from |53j . For 
the WMAP cosmological parameters 3 , we obtained the 
lower bound 



n > 4.2 X 10^ (95% C.L.). 



(15) 



As we expected from Fig. [U the kinetic gravity braid- 
ing model with small n is obviously ruled out, while the 



^ Most ISW detections have less than 3(7 detection level |49|473|| . 
In addition, the authors in Ref. [zsj claimed the possibility of 
the absence of ISW signal extracted from the five-year WMAP 
data and the AAOmega LRG redshift survey. Although the de- 
tection of the cross-correlation between the cosmic microwave 
background and the large scale structure has been under debate 
(see, e.g., ff^). we use the result of a combined analysis of the six 
catalogues in Ref. [s^l . which is the highest detection (~ 4.5cr). 



model with large n is favored by observations. This is 
because the ISW effect anti-correlates with the LSS for 
n < O(IOOO); and the sign of the cross-correlation func- 
tion determines the lower bound (jTS]). Note that the 
galaxy bias in Ref. [s^ is obtained from the two-point 
correlation functions of each catalog assuming the ACDM 
model. In our case, the matter power spectrum for 
n > 1000 and 0.001 /iMpc"^ < fc < 0.1 hMpc~^, which 



is the dominant wavenumbcr in Eq. (jl3p . is almost the 
same as those in the ACDM model. Therefore, this is not 
problematic, and we also confirmed that using the galaxy 
bias of the ACDM model does not significantly change 
the results of chi-squared analysis. The lower bound ([15]) 
changes within only 10 percent even when we assume the 
bias as free parameters. This is because the constraint on 
the model parameter n is determined by the sign of the 
cross-correlation function, which is almost independent 
of nature of bias. 



V. PARAMETRIZED MODEL 

In the previous section, we have demonstrated that 
the ISW-galaxy cross-correlation strongly constrains ki- 
netic gravity braiding. The key to this result was the 
time-evolution of the gravitational potential controlled 
by the effective gravitational coupling Gofr in the den- 
sity perturbation equation. In this section, we employ a 
phenomenological paramctrization of the time- evolution 
of Goft rather than Goff derived from a concrete model, 
and investigate how the parameters are constrained by 
tomographic observations. This approach is useful to es- 
timate a general prospect of constraints on generalized 
models (e.g., [36j). 

Provided that the quasi-static approximation is valid, 
the evolution of overdensities is described by the equation 
of the form ^ with ^ in all the scalar-tensor theories 
with second-order field equations [36| including kinetic 
gravity braiding. We parametrize the evolution of Goff 
as follows: 

Goff 



G 



1 



92 



(16) 



where gi and (72 are the free parameters. This repro- 
duces the effective gravitational coupling e.g., of the ki- 
netic gravity braiding model considered in this paper for 
small n. In particular, the effective gravitational coupling 
with gi = 1 and 32 = 4 yields that of the kinetic gravity 
braiding model for n = 1 within a few percent precision. 
We consider the model whose sound speed for the scalar 
field does not approach zero. This means the quasi-static 
approximation is applicable for the dominant wavenum- 
bers in Eq. Then, the evolution equation for the 

matter overdensity and the Poisson equation are of the 
form Eqs. ([5]) and (0). In the case of a minimally cou- 
pled scalar field model, the traceless part of the Einstein 
equations gives 4* -I- $ = 0. 

In this section, wc present the observational constraints 
on the above parametrized model and perform the Fisher 
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FIG. 1: The cross-correlation function theoretically calculated by using Eq. (|12p and the data obtained in [53]. Each curve shows 
the cross-correlation function of the ACDM model (solid curve), the KGB model (dashed curve) with n = 5000, 1000, 100, 10, 
and 1, from the top to the bottom, respectively. 



matrix analysis. Figure [2] shows the contour of Ax^ on 
the gi — 52 plane. We assumed that the background 
expansion history is the same as that of the standard 
ACDM model, and used the cosmological parameters 
from the WMAP7 data 0. The background history 
for n > 100 is almost identical to the LCDM model. 
As seen from the section IV, for n < 100 the ISW-LSS 
cross-correlation is negative, which is clearly incompati- 
ble with observational data, and the different background 
evolution does not change the sign of the correlation 
function and hence is not important. For each param- 
eter, we assume that the galaxy bias is determined by 
the amplitude of the power spectrum of the matter dis- 
tribution, = 6(ACDM)^^(4«))/^(acdm)^^(0^^ ^^^^^ 

^(ACDM) ^Yie bias used in the previous section, /j^^'-'^^) 

is the growth factor of the ACDM model, and z**-* is 
the mean redshift of the i-th catalog. Our result implies 
that the deviation of the effective gravitational coupling 
from Newton's constant must be very small, as expected 
in [t^]. Since gi — > or 32 — >■ oo corresponds to the 
ACDM model for a < I the constraint on f/2 can not be 
determined as long as the ACDM model is favored by 
observations. 

Next, we give a forecast for the accuracy of constrain- 
ing gi and 172 from future galaxy surveys. In the Fisher 
matrix analysis [78| . we consider the following redshift 



distribution of the galaxy sample per unit solid angle, 

/51 



dN _ Ngl3 
d^"zo"+ir((a + l)//3) 



z exp 



where a, f3, zq are the parameters, and Ng 
The mean redshift is determined by 



1 



dzz 



dN zoT{{a + 2)/ 13) 



dz 



T{{a + l)/f3) 



(17) 



/ dN/dz. 



(18) 



We used a = 0.5, /3 = 3, iVg = 35 arcmin^^ and z„i = 
0.9. We divide the galaxy sample into four redshift bins, 
0.05 < z < 0.6zm, 0.6z„i < z < Z„i, z„i < z < \AZm, 
and lAzm < z < 2.5 [t^, where the redshift distributions 
of i-th bin (zi-i < z < Zi) arc given by [soj 



W,(z) 



Ai dN 
Y~dz 



erfc 



V V2a{z) 



erfc 



\V2a{z) 



(19) 



where erfc is the complementary error function, Ai is de- 
termined by the normalization constant, and we adopted 
cr(z) = 0.03(1 -f z). The Fisher matrix is given by [811483} 



i^^^=/sky^(2/+l) 
I 



56, 



,(20) 



5 
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-0.4 




2 4 6 8 10 
gi 

FIG. 2: Contour of Ax^ on the g\ — g2 plane. The dashed 
curve and sohd curve are the la and 2a confidence contour- 
levels, respectively. 

where /sky is the fraction of sky common to both the 
CMB and the galaxy survey maps, (^^^w-G -g ^j^^ 
point angular cross-correlation defined by Eq. 8^ is 
the model parameter gi and g2, and 

COV(0 = (Cf W-G)2 ^ (^ISW ^ ArCMB)(^G ^ 7Vg).(21) 

Figure [3] shows the 1-sigma (dashed curve) and 2-sigma 
(solid curve) confidence contours in the gi — (?2 plane, 
assuming future survey such as the HSC-like weak lens- 
ing survey {A = 1500 deg^) and the LSST-like survey 
{A = 20000 deg'^). The target model is the parametrized 
model given by Eq. (|T6l) with gi ~ 0.1 and 52 = 3. Wc 
used the WMAP cosmological parameters for iloh^, ^bh^ 
and Atj and assumed that the bias and the initial am- 
plitude of the fiuctuation is determined by individual ob- 
servations such as CMB power spectra and matter power 
spectra. Then we consider only 2 parameters gi and 
(72 in the Fisher matrix analysis. As one can see from 
Fig. 131 the tomography, by dividing the galaxy sample 
into several redshift bins, can significantly improve con- 
straints on the parameter 32- Therefore, the ISW-LSS 
cross-correlation with tomography can dramatically im- 
prove the determination of the evolution of the effective 
gravitational coupling as well as its amplitude. 

VI. CONCLUSION 

In this paper, we have obtained observational con- 
strains on the kinetic gravity braiding model from the 
cross-correlation function between the ISW effect and 
the large scale structure. We have found that the cross- 
correlation function of the kinetic gravity braiding model 
with small n is negative, and therefore such a model is 
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0.2 /■ /' 

/ ' _ _ — ✓ / 




-0.1 : 

2 4 6 8 10 
g2 

FIG. 3: The 1-sigma and 2-sigma contours in the g\ — g2 
plane for the future survey. We assumed A — 1500 deg^ + 
4 bins (dot-dashed curve), A — 20000 deg^ (dashed curve), 
A = 20000 deg^ 4- 4 bins (solid curve). The target parameter 
is gi =0.1 and g2 = 3. 



not favored by observations. Constraints from ISW-LSS 
cross-correlation is stronger than those obtained from 
SNIa and the CMB shift parameter Thus, the co- 

variant galilcon model, n = 1, is ruled out by the cross- 
correlation between the ISW effect and the large scale 
structure. 

In deriving the constraints, we have clarified that the 
evolution of the effective gravitational coupling is crucial 
for the behavior of the ISW effect. Motivated by this 
fact, we also considered the model with a phenomeno- 
logical parametrization of Gcff that can reproduce the 
effective gravitational coupling of the covariant galileon 
model. We have confirmed that the deviation from the 
cosmological constant must be very small: the 1-sigma 
(2-sigma) error in determining 171 is Agi ^ 0.15 (0.3) 
when g2 ^ 2. We have performed the Fisher matrix anal- 
ysis, assuming a HSC-like survey and a LSST-like survey, 
and found that the ISW-LSS cross-correlation with the 
sample divided into multiple redshift bins can improve 
significantly the constraint on g2- This tells us that the 
ISW-LSS cross-correlation obtained from the future sur- 
vey would be a powerful probe for long distance modifi- 
cation of general relativity. 
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